Comparison between the ages (after controlling for type and gender)

Questions

  • What are the differences between the ages?
  • Which genes and pathways are differentially expressed between 8w and 52w, between 52w and 104w, between 8w and 104w?
  • Are they the same? Is there a gradient?

Age effect - General Questions

  • What are the differences between the ages?
  • Which genes and pathways are differentially expressed between 8w and 52w, between 52w and 104w, between 8w and 104w? Are they the same? Is there a gradient?
  • Are they different for the two genders?
  • Are they different for the two types?

Loads

Libraries and functions

Bioconductor version 3.6 (BiocInstaller 1.28.0), ?biocLite for help

==========================================================================
*
*  Package WGCNA 1.63 loaded.
*
*    Important note: It appears that your system supports multi-threading,
*    but it is not enabled within WGCNA in R. 
*    To allow multi-threading within WGCNA with all available cores, use 
*
*          allowWGCNAThreads()
*
*    within R. Use disableWGCNAThreads() to disable threading if necessary.
*    Alternatively, set the following environment variable on your system:
*
*          ALLOW_WGCNA_THREADS=<number_of_processors>
*
*    for example 
*
*          ALLOW_WGCNA_THREADS=4
*
*    To set the environment variable in linux bash shell, type 
*
*           export ALLOW_WGCNA_THREADS=4
*
*     before running R. Other operating systems or shells will
*     have a similar command to achieve the same aim.
*
==========================================================================


Allowing multi-threading with up to 4 threads.

preparing gene to GO mapping data...
preparing IC data...
preparing gene to GO mapping data...
preparing IC data...
preparing gene to GO mapping data...
preparing IC data...

Data

Stats

Wald padj < 0.05LFC > 0 (Wald padj < 0.05)LFC < 0 (Wald padj < 0.05)
52w VS 8w21571176 981
104w VS 8w319217011491
104w VS 52w228412521032

Differentially expressed genes

Warning message in grep(txdbPattern, installedPackages):
“argument 'pattern' has length > 1 and only the first element will be used”Warning message in pcls(G):
“initial point very close to some inequality constraints”Warning message in grep(txdbPattern, installedPackages):
“argument 'pattern' has length > 1 and only the first element will be used”Warning message in pcls(G):
“initial point very close to some inequality constraints”Warning message in grep(txdbPattern, installedPackages):
“argument 'pattern' has length > 1 and only the first element will be used”Warning message in pcls(G):
“initial point very close to some inequality constraints”

Counts

Genes that are differentially expressed in 52W vs 8W with abs(FC) > 2

Genes that are differentially expressed in 104W vs 8W with abs(FC) > 2

Genes that are differentially expressed in 104W vs 52W with abs(FC) > 2

Compare the numbers

Differentially expressed genes

Differentially more expressed genes

Differentially less expressed genes

Compare the intersections in term of expressions

Genes that are differentially expressed in 52W vs 8W and in 104W vs 8W (and not in 104w vs 52w)

52w > 8w & 104w > 8w52w > 8w & 104w < 8w52w < 8w & 104w > 8w52w < 8w & 104w < 8w
5630 0 329

With abs(FC) > 1

Genes that are differentially expressed in 52W vs 8W and in 104W vs 52W (and not in 104w vs 52w)

52w > 8w & 104w > 52w52w > 8w & 104w < 52w52w < 8w & 104w > 52w52w < 8w & 104w < 52w
0 1072750

With abs(FC) > 1

Genes that are differentially expressed in 104W vs 8W and in 104W vs 52W (and not in 52w vs 8w)

104w > 8w & 104w > 52w104w > 8w & 104w < 52w104w < 8w & 104w > 52w104w < 8w & 104w < 52w
5040 0 567

With abs(FC) > 1

Genes that are differentially expressed in 52W vs 8W, in 104W vs 8W and in 104W vs 52W

52w > 8w & 104w > 8w & 104w > 52w52w > 8w & 104w > 8w & 104w < 52w52w > 8w & 104w < 8w & 104w > 52w52w > 8w & 104w < 8w & 104w < 52w52w < 8w & 104w > 8w & 104w > 52w52w < 8w & 104w > 8w & 104w < 52w52w < 8w & 104w < 8w & 104w > 52w52w < 8w & 104w < 8w & 104w < 52w
1896 0 7 15 0 6 46

With abs(FC) > 1

DEG into gene co-expression network

  • White: up-regulated
  • Black: down-regulated

Order: 52w VS 8w, 104w VS 52w, 104w VS 8w

GO analysis

Biological process

Dot-plot with the most over-represented GO (20 most significant p-values for the different comparison)

Network based on description similarity

52w VS 8w

scatterplotThree

104w VS 52w

scatterplotThree

104w VS 8w

scatterplotThree

Cellular components

Dot-plot with the most over-represented GO (20 most significant p-values for the different comparison)

Molecular functions

Dot-plot with the most over-represented GO (20 most significant p-values for the different comparison)

KEGG pathways

[1] "Note: 7 of 4750 unique input IDs unmapped."
[1] "Note: 7 of 4750 unique input IDs unmapped."
[1] "Note: 7 of 4750 unique input IDs unmapped."
[1] "Note: 7 of 4750 unique input IDs unmapped."
[1] "Note: 7 of 4750 unique input IDs unmapped."
[1] "Note: 7 of 4750 unique input IDs unmapped."
[1] "Note: 7 of 4750 unique input IDs unmapped."
[1] "Note: 7 of 4750 unique input IDs unmapped."
[1] "Note: 7 of 4750 unique input IDs unmapped."
[1] "Note: 7 of 4750 unique input IDs unmapped."
[1] "Note: 7 of 4750 unique input IDs unmapped."
[1] "Note: 7 of 4750 unique input IDs unmapped."
[1] "Note: 7 of 4750 unique input IDs unmapped."
[1] "Note: 7 of 4750 unique input IDs unmapped."

Pathway graphs available at ../results/dge/age-effect/age/over_repr_kegg/

[1] "Note: 7 of 4750 unique input IDs unmapped."
[1] "Note: 7 of 4750 unique input IDs unmapped."
[1] "Note: 7 of 4750 unique input IDs unmapped."

Warning message in structure(x$children, class = "XMLNodeList"):
“Calling 'structure(NULL, *)' is deprecated, as NULL cannot have attributes.
  Consider 'structure(list(), *)' instead.”Warning message in structure(x$children, class = "XMLNodeList"):
“Calling 'structure(NULL, *)' is deprecated, as NULL cannot have attributes.
  Consider 'structure(list(), *)' instead.”Warning message in structure(x$children, class = "XMLNodeList"):
“Calling 'structure(NULL, *)' is deprecated, as NULL cannot have attributes.
  Consider 'structure(list(), *)' instead.”Warning message in structure(x$children, class = "XMLNodeList"):
“Calling 'structure(NULL, *)' is deprecated, as NULL cannot have attributes.
  Consider 'structure(list(), *)' instead.”Warning message in structure(x$children, class = "XMLNodeList"):
“Calling 'structure(NULL, *)' is deprecated, as NULL cannot have attributes.
  Consider 'structure(list(), *)' instead.”Warning message in structure(x$children, class = "XMLNodeList"):
“Calling 'structure(NULL, *)' is deprecated, as NULL cannot have attributes.
  Consider 'structure(list(), *)' instead.”Warning message in structure(x$children, class = "XMLNodeList"):
“Calling 'structure(NULL, *)' is deprecated, as NULL cannot have attributes.
  Consider 'structure(list(), *)' instead.”Warning message in structure(x$children, class = "XMLNodeList"):
“Calling 'structure(NULL, *)' is deprecated, as NULL cannot have attributes.
  Consider 'structure(list(), *)' instead.”

Pathway graphs available at ../results/dge/age-effect/age/under_repr_kegg/